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I. OVERVIEW 


This report is basically a compendium of memos submitted to Dr. A. 
Chang and Dr. D. LeVine of NASA/GSFC during the course of this contract. Also 
included are copies of the vugraphs presented by J. Carr of ORI to the NASA 
Earth Observing System/High resolution Multifrequency Microwave Radiometer 
(EOS/HMMR) panel. This presentation covers, our preliminary work with aperture 
synthesis concepts for EOS, and appeared in our preliminary report. New to 
the final report is an exposition of the effects of nonvanishing bandwidths on 
image reconstruction in apperture synthesis systems. It is found that 
nonvanishing bandwidths introduce errors in off-axis pixels when naive Fourier 
processing is used. The net effect is for bandwidth to limit sensor 
f ield-of-view. Although known to radio astronomers, this effect has to our 
knowledge never previously been quantified. To quantify this effect we wrote 
a computer program (a copy of which is delivered with this report) which is 
documented in this report. Example runs are included which illustrate the 
resultant radiometric errors and effective f ields-of-view for a plausible 
simple sensor. 



II. REAL APERTURE SYSTEMS 


The study began by considering real aperture systems for EOS. These 
concepts grew out of the considerations for a Shuttle-based instrument 
examined under a previous contract with ORI. 

It was decided that the instrument concepts should address the 
measurement requirements of: 

1) 10 km spatial resolution 

2) 3 day temporal resolution 

3) 1 K radiometric accuracy. 

The requirements were derived from the "EOS Science and Mission Requirements 
Working Group Report" (NASA-TM-86129) and the published experience of other 
instruments (e.g., Skylab S-194). 

The first memo attached considers an electronically scanned array 
with a single beam. It was found that the constraints on integration time by 
the need to form an image did not permit all of the measurement requirements 
to be simultaneously satisfied. A factor of two improvement in the 
sensitivity is required to meet .5 K radiometric sensitivity. The remainder 
of the 1 K radiometric accuracy is budgeted for systematic errors. Several 
ways to enhance performance were then studied. These are: 

1) Increase the reception bandwidth beyond the 27 MHz reserved for 
passive use only 

2) Build a total power radiometer 
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3) Build a reference averaging radiometer 

4) Build a multi beam (wi skbroom) radiometer. 

The second attached memo examines the feasibility cf baniwidth 
extension. It is rejected due to the large number and nigh powers of 
interferers in neighboring bands. The third memo deals with total power 
radiometry. It appears that total power radiometry may be a feasible option, 
but that temperature and power supply stability are crucial. On the 
performance spectrum between total power and Dicke radiometers is the 
reference averaging radiometer. This is described in the fourth memo and does 
appear attractive. 

Finally, a wiskbroom system with four beams will certainly provide 
the necessary sensitivity; however, this introduces antenna design 
complications. Four beams are necessary since the number of beams increases 
the integration time linearly and the sensitivity is inversely proportional to 
the square-root of the integration time; hence the factor of four is required 
to make-up the factor of two shortfall in sensitivity. 
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January 25, 1985 


TO: File 

FROM: 0. Carr and B. Candey 

RE: Bandwidth Extension 

Analysis to establish a baseline real aperture L-band sensor for EOS 
has found that Integration time constraints do not permit both the desired 
spatial resolution and the desired radiometric resolution. It Is also 
anticipated that this Integration time constraint will place more severe 
restrictions upon synthetic aperture systems. One possible solution is to 
Increase the bandwidth from the nominal 27 MHz (1400-1427 MHz). Initially 
this seemed feasable since the 1370-1400 MHz band Is designated for passive 
use on a secondary basis. However, as the following analysis will indicate 
bandwidth extension Is infeasable due to the presence of numerous high powered 
transmitters. 

The frequency bands neighboring 1400 MHz have the following 
allocations^: 

1) The 1300-1350 MHz band is allocated to radionavigation worldwide 
except for fixed and mobile services in Indonesia and radiolocation 
service In Ireland and the United Kingdom. Radionavigation use is 
limited to ground-based radars and associated airborne transponders. 
The U.S. government allocations also allow radiolocation but limited 

l 

to the military services. 


^ Manual of Regulations and Procedures for Federal Radio Frequency 
Management , NTIA, May 1984. p. 4-53ff. 



2) The 1350-1400 MHz band Is allocated to fixed, mobile, and 
radiolocation services In Region 1 (North and South America) and 
radiolocation In Regions 2 and 3 (Europe, Africa and Asia) with some 
existing radionavigation installations In the Eastern Bloc, U.S., and 
Canada. The 1370-1400 MHz band is also allocated to passive space 
and earth research and observations on a secondary basis. In 
addition, the frequency 1381 MHz Is allocated In the US to Fixed and 
Mobile Satellite Services (Space-to-Earth) for the relay of nuclear 
burst data. In the U.S., over 100 high powered radars operate In the 
1370-1400 MHz band and International listings also Indicate 
extensive use. 

3) The 1400-1427 MHz band Is reserved solely for passive radio 
astronomy, space research, and earth exploration-satellite 
observations. Some second and third harmonics and various 
Intermodulation products of television transmitters and extraband 
radiation from radars may fall within the 1400-1427 MHz band. 

4) The 1427-1525 MHz bands are allocated to fixed, mobile, and 
space operation and are used heavily. 

Ideally one wants a large bandwidth for enhanced sensitivity; 
however, artificial sources are many orders of magnitude brighter than the 
thermal emmisslons of the earth and must be excluded from the Instrument 
bandpass. We determine the interference threshold as follows. The required 
receiver sensitivity is &T - .5 K which corresponds to a power level of k-OfifO'B 
with It = 1,38 x IQ" 23 J/K and B the bandwidth. The maximum sensor- 
Interferer range will occur just when the sensor rises above the Interferer's 
horizon (see Figure 1). Coincidently this Is where radiolocation services 

t 

(e.g., FAA radars) would have their maximum gain directions and so this Is the 
anticipated worst case geometry dispite the large range. At a nominal 
altitude of h = 700 Km, the sensor - interferer range is R = 3069 Km with the 
Interferer at a 64° angle from the sensor's nadir. We assume a nominal 10m 
x 10m sensor aperture which scans A 33° to cover a 900 Km swath. This 
places the interferer 31° off the maximum gain direction at the extreme 
point of the scan. A uniformly weighted aperture would have a power pattern 
which falls as 

(x -‘♦y* 

l4 
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Figure 1. Sensor - Interferer Geometry 


where L Is the aperture dimension (10m), X Is the wavele- gth (21 cm), and 
Is the angle off the maximum gain direction. For $ * 31° thli puts the peak 
sldelobe level at 6 * -38 dB. However, It Is better to taper the aperature 
weighting In order to trade beam width for lower sldelobes (1r creased main 
beam efficiency). Sldelobe levels at 31° for a 90X main beam efficient 
antenna would be more like G * -50 dB (see Figure 2). The power Into the 
sensor's receiver from the Interferer Is then 

p 

r — . 0) »- ^ 


where EIRP is the Interferer's Equivalent Isotropic Radiated Power. In which 
case, the maximum tolerable EIRP to not exceed & power threshold P Is given by 


(Elfc?) = 


4ttK. z P 


We will conservatively set P = k(£T)B which Is the receiver 
sensitivity. For B = 27 MHz, P = -157 dBW and for B * 100 MHz, P = -152 dBW. 
Using R * 3069 Km, L = 10 m, G c -50 dB gives EIRP thresholds of + 14 dBW and 
+19 dBW for the 27 MHz and 100 MHz bandwidths respectively. 


Typical EIRPs for fixed stations operating near 1.4 GHz are 37 dBW to 

o 

55 dBW c which would definitely represent an interference threat. 

Additional sidelobe rejection of -36 dB would be required to place the 55 dBW 
transmitters at the received power threshold -152 dBW, and another -10 dB 
should be added to place the Interferer below the recleved power threshold 
with adequate margin. Then a sidelobe level of -96 dB is required at 31° 
off the maximum gain direction. This would be extremely difficult to achieve. 


^Recommendations and Reports of the CCIR, 1982: XVth Plenary Assembly in 

6eneva, Volume II: Space Research and Radioastronomy, ITU, p. 386. 
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Figure 2. Power Pattern for a 90% Efficient Antenna. 

(Source: "L-band Radiometer Antenna Study", Hughes Aircraft 

under contract NAS5-23418, May 1977) 





Using a received power threshold of -165 dBK (compare with -157 dBW) 
corresponding to an EIRP threshold of 12 dBW (compare with 19 dBW), the 
International Radio Consultative Committee (CCIR) has performed an 
Interference analysis for a 100 MHz bandwidth sensor 3 . The results are 
shown In Figure 3 where total data loss occurs In the shaded regions. This 
analysis Includes only the fixed services which populate the 1350-1400 MHz 
band, so the actual Interference environment Is In fact worse. 

Transmitters In the 1250-1400 MHz band Include FAA aeronautical 
radionavigation radars, DoD long range surveillance and search radars, DoD 
communication radios, the GPS downlink for relay of nuclear burst data, the 
TLPS Army transponder and classified military transmitters. A conmon FAA 
radar, the ARSR-1, operates at frequencies of 1280-1350 MHz, has an antenna 
gain of 34.2 dBi and 4 MW peak power over 2 microsecond pulses. DoD 
communication radios, such as the AN/GRC-50 and AN/GRC-103, broadcast In the 
1350-1849.5 MHz band with antenna gain of up to 20 dBi and peak power of 8-30 
W. The TLPS transponder uses the 1350-1400 MHz band with a peak power of 120 
W.« 


In light of the above interference analyses, it appears that 
bandwidth extension beyond the reserved 1400-1427 MHz band Is infeasable. 

This conclusion is corrobated by the opinion of David P. Struba, NASA 
Headquarters Frequency Management Program Manager, that the secondary 
allocation for passive use of the 1370-1400 MHz band Is almost 
meaningless. ^ In addition there Is also concern that extraband radiation 
and TV transmitter harmonics could pollute the 1400-1427 MHz band. Although 
there is no mention of interference problems in the literature reporting the 
1973 Skylab L-band radiometer results, this Issue should be studied further. 

If it is determined that such emissions might pose an Interference threat then 

i 

a Shuttle survey with a low cost low resolution antenna might be considered to 
characterize the L-band environment. 

3 Ibid . pp. 379 - 399. 

^Andrew Farrar, Spectrum Resource Assessment In the 1215-1400 MHz Band , NTIA 
Report 81-83, Sept. 1981, pp. 36-8, 50. 

^Private communications between David P. Struba and Bob Candey. 
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TO: 
FROM: 
SUBJ : 


J. Carr 

R. E. Prince f&P 

L - Band Amplifier Drifts with Respect to Total 


Power Radiometers 


I have made a number of phone calls and managed to talk to about half 
of the people I attempted to contact. 

Two attempted contacts that never returned their calls were: 

Chase Hearne, a specialist In GaAs amplifiers. 

Hans Blume, a specialist In uncooled parameter amplifiers. 

They are both NASA/Langley (LaRC) people whose names were given me by 
their boss at LaRC, Richard Harrington (804) 865-3631. Harrington said they 
would call, but ... So be It. I mention this only because these two were 
supposedly hot numbers (technically speaking). I got Harrington's name from 
John Shlen at GSFC. 


A man named Joe Rogers (GSFC, 344-8816) has tested a number of GaAs L 
r - band amps. He particularly likes units manufactured by California Amplifier 
and by Berkshire Technologies, Inc. (also in California). Rogers indicates 
that another GSFC person, David Buhl (344-8810) has a noise figure and gain 
test setup for these devices based on an HP 8970 Instrument. Whether he can 
measure the small variations we are looking for Is another question. Buhl 
also suggests that we talk to Robert Jones, RF Technology Branch in GSFC, Bldg 
19. Apparently building 19 has a test setup based on an AIL Instruments. 


Although nobody I talked to had test or analytical data for small 
scale, short term 6aAs L - band amplifier gain stability with temperature, all 
my correspondents agreed that primarily temperature (and only slightly supply 
voltage) controls amplifier gain. All agreed that with sufficient temperature 
stabilization a gain stability of 10"3 (-30 db) was quite feasible. 

Supply voltage requirements are typically +_ 15 volts § 50 ma. 
Regulation should be # 10% or better. 


All of the amplifier people I spoke to quoted specifications for an 
L - band unit similar to the following: 

frequency range: 1.30 to 1.60 GHz 

noise figure: .7 to 1.0 db at room temperature 

aNF/^T approx. .05 db/10°c (down to a noise figure of 
approx. .6 at -70°c) 

gain: 35 to 40 db (power) at room temperature 
aG/^T: -.02 db/Oc to -.04 db/°c (gain increases with 

dropping temperature) 
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This /sG/^T number corresponds to approximately .1° C temperature 
stability to meet a -30 db *sG/G (or .01° C stability for a -40 db aG/G). 
0See attached calculation). 

I have enclosed some test data I received from California Amplifier, 

Inc. 


I talked to Jim Bremer for a short while on the telephone concerning 
his reference - averaging radiometer design. He Indicates that In 
approximately 1977, Hughes Alrcarft did a computer simulation and then built a 
labatory breadboard of the reference - averaging technique. Jim said that the 
Hughes work verified the Bremer calculations but that simultaneous advances in 
amplifier stability mode a total-power system feasible. The Hughes work was 
at millimeter wave lengths. 

Jim brought up something else to be remembered about total power 
systems: a D.C. - Restore capability probably will be necessary. D.C. - 

Restore Is essentially is a calibration - it returns the analog output to the 
A/D converter to a known D.C. level for a given input every now and then. 
Typically, It is done once per scan line (every few seconds) or whenever 
practical. If the amplifier Is stable enough that it's baseline plus dynamic 
range never clamps or extends beyond the A/D converter useful range of Inputs, 
the D.C. - Restore function as a specific analog circuit could be omitted. 

This function could then be subsumed under the calibration function and 
accomplished digitally at the output of the A/D converter (or still as an 
analog function, of course). 



People I have attempted to contact Include: 


person and Location Successful 


R. Harrington, LaRC x 

(804) 865-3631 
Chase Hearne, LaRC 
Hans Blume, LaRC 

31m Shleu, GSFC x 

Roger Ratliff, GSFC x 

Joe Rogers, GSFC x 

David Buhl, GSFC x 

Trontech Corp., Neptune, N.J. 

(201) 229-4348 

Miteq, Hauppage, N.Y. x 

(Jared Slddiqul) (516) 543-8873 
Applied Microwave Corp., Lawrence, Kansas x 

(Dr. David Brunfleld, formerly U.KS Remote 
Sensing Center) (913) 749-3511 
California Amplifier, Newbury Park, Cal. x 

(John Ramsey) (805) 499-8535 
Berkshire Technologies, Inc., Oakland, Cal. x 

(415) 655-1986 


Unsuccessful 


x 

x 


x 


Although I did not call them, several people suggested Amplica, 
.(805) 498-9671. They are a division of COMSAT. 
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III. APERTURE SYNTHESIS SYSTEMS 


As an alternative to real aperture systems we have considered 
synthetic aperture systems. Aperture synthesis can offer high spatial 
resolution with a sparse antenna array and so it can be economically 
attractive. Initially the more general case of a synthesis/scanning hybrid 
was considered (of which a pure synthesis instrument is a degenerate case). 
Such a hybrid would step-scan its f ield-of-view across-track. Aperture 
synthesis provides the imaging within each field-of-view. The first of the 
following memos describes a one-dimensional scanning/synthesis hybrid. It 
concludes that such instruments can perform the soil moisture mission. An 
especially remarkable result is that the necessary radiometric resolution 
could be achieved despite the diminished aperture area. This is because the 
collection area loss is offset by both an increased integration time (set by 
the real f ield-of-view) and the parallel channel structure of an aperture 
synthesis instrument. 

The second of the following memos details a two-dimensional 
synthesis/scanning hybrid concept. It concludes that the sensitivity of such 
an instrument is identical to the one-dimensional case for given aperture 
area. Therefore, there is no advantage to this more complicated scheme and 
only one-dimensional concepts were henceforth considered. 

Following this memo is a set of notes detailing the theory of 
aperture synthesis. It covers the signal processing theory, resolution, 
grating lobes, and nonlinear effects (fringe washing). The grating lobes and 
fringe washing effects are particularly important as they restrict the field- 
of-view. 


if 



Fin< 11y, details of how to construct a pure aperture synthesis 
instrument we^e considered. In the last three memos an analog and two digital 
architectures are [resented. All architectures have the same external data 
rates, wherea; the digital architectures (they do correl lations digitally as 
opposed to using analog mixers) can have very significant internal rates. 
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IV. EOS/HMMR PANEL PRESENTATION 


What follows are the vugraphs from our presentation to the EOS/HMMR 
panel. This presentation summarizes some of our results and was intended to 
generate some confidence that aperture synthesis could be used to perform the 
soil moisture mission. We took as science requirements: 

1) 10 km spatial resolution 

2) 3 day temporal resolution 

3) 1 K radiometric resolution. 

We took the 3 day temporal resolution requirement to mean that both 
dayside and nightside overpasses count as revisits. At a presumed 700 Km 
sun-synchronous circular orbit, a 450 km swath implies three day coverage (at 
least in a statistical sense, i.e., mean revisit time). Since swath width is 
the most difficult parameter to achieve due to grating lobes and fringe 
washing, and since we wanted to present a "minimal" instrument, the swath 
width goal was degraded to 400 km. This implies 3 day (statistical) coverage 
in the temperate latitudes above 30° and 4 day (statistical) coverage at the 
equator. Increasing the swath width to 450 km (or even 900 km which gives 
precisely periodic 3 day coverage rather than statistical 3 day coverage) is 
possible at the cost of more antennas. 

The system presented consists of 11 13.2 m x .34 m antennas occupying 
a 13.2 m x 16.2 m array. The array fill-factor is 23 percent. It provides 10 
km resolution both along and across-track with a .6 K radiometric 
sensitivity. With this one-dimensional concept, aperture synthesis provides 
the across-track resolution; whereas, the width along-track of the antennas' 
fields-of-view provides the along-track resolution. 



A digital system architecture was presented. This arcnitecture joins 
11 antenna/receiver units, arranged to form a minimum redundancy linear array, 
with about 40 in-phase and 40 quadrature digital correlators. Each pair of 
an ennas forming a unique baseline are joined with both in-phase and 
quadrature correlators. In addition, several autocorrelation channels would 
be required to measure the zero-baseline or DC component of the scene spatial 
frequency distribution. Each channel need not physically correspond to a 
unique piece of hardware, as different baselines could be time-phase 
multiplexed on a common correlator. The instrument internal data rate would 
be fairly high at 20 Mbps which is determined by the Nyquist sampling rate of 
the signals being digitized by the antenna/receiver units. The external data 
rates would, however, be quite modest at about 1 Kbps which is determined by 
the dynamic range and noise levels in the measurements. Fourier processing to 
construct an image could be easily done on the ground in real-time. 
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TEMPORAL RESOLUTE AND SWATH WIDTH 
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SOIL HOISTWRE-BRIGHTJESjlt«PERATllRE RELATIONSHIP 
(As Reported for Sky lab 8-19*0 
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APERTURE SYNTHESIS 
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Individual Antenna Patterns Determine Field-of-View 
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Integration Time Is Set As Time Required To Advance One 
Resolution Element Along-Track 
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Necessary Sensitivity Achieved With Fill-Factor 
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Radiometric Resolution: 1 K 
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0RK3HAS, !PA@E SB 

OF POOR QUAlffiY 


V. BANDWIDTH AND RECONSTRUCTION 


As noted earlier, a nonvanishing system bandwidth results in imaging 
errors ("fringe washing") off-axis when naive Fourier processing is used for 
image reconstruction. The net effect is for system bandwidth to limit 
f ield-of-view. In this section we exhibit the theory behind this effect. 

This is exploited in the computer simulation described in section VI. 

The first attached memo derives the relationship between system 
response and scene brightness distribution for a polychromatic extended 
source. The problem of image reconstruction is cast as inverting an integral 
equation: 

fOO 

r r /i ziri ux 

= J J kco>)\<(x)£ ° dxcJeo 

o -oo 

where 

R(u) is the response at spatial frequency (or baseline) u 

G(u>) is the filter power gain at radian frequency^) 

u ) 0 is the center or fiducial radian frequency 
K(x) is the integral kernal to be estimated at scene coordinate x 
(sine of angle across track). 
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ORIGIN PASE SS 
OF POOR QUALIFY 

The kernal K(x) is the product of the antenna power pattern, scene bnghtne.s 
distribution and Jacobian factor (i—X 2 ) ^ This is almost a Fourier 
relationship between R(u'i and K(x). 


To see the connection between bandwidth and f ield-of-view, we perfo~m 
a change of variables = (ca-co„ ) and integrate over ^ . Then 

<-°o 

+cO 

F , ZTlK-X , 

=• J ku.,x) «Tx)£ dx 

- va 


where 


- 



-riO 

j co 0 G[('-+x)^c) 5 > -^ 

jo ^ - 1 


Here h(u,x) multiplies K(x) which contains the antenna power pattern. Then 
h(u,x) can be thought of as an additional but spatial frequency dependent 
pattern factor. For a single frequency bandpass h(u,x) in a constant, but as 
the bandwidth increases, it becomes tapered away from x=0. 


The second memo describes how this formalism is translated into the 
discrete problem of inverting a finite sampling of spatial frequencies to 
estimate scene brightness. There we have an almost discrete Fourier 
relationship: 

2TT t (_m — i )(.v\ - i V ^ 

Z k„vj nM e 

n = i 


where 


R m is R ( u) sampled at and weighted by a phase factor 

K p is K(x) sampled at x = ( n - x 

W nm is a filter design dependent weighting 


lil 



The samplings A\, Au- are such that A\ Aia - ^ where N s the 


number of baselines. 


For a real sensor some baselines may occur with multiplicities other 
than 1, in which case, we multiply by the multiplicity factor 
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VI. COMPUTER SIMULATION 


This section documents our aperture synthesis computer simulation 
designed to quantify the nonvanishing bandwidth and reconstruction problem. 
The source listing follows, it is written in Microsoft FORTRAN for the IBM 
PC. The source and object codes for which are delivered on a floppy disk. 

Two main programs appear - APSYNSIM and BASELINE. APSYNSIM calls subprograms 
FOURT, SINC, and PLOT. 

APSYNSIM is the main program for the aperture synthesis simulation. 
It reads the run input data from file MASTER.DAT (listed after source 
listings). The inputs are as in table VI-1. System parameters specified are 
the number of baselines, null-to-null field-of-view, center frequency, and 
bandwidth. The test scene is assumed to be of constant nominal brightness 
temperature and limited to a specified scene width (<180°) to avoid the 
Jacobian factor in K(x) raising a division by zero error. APSYNSIM reads the 
baseline multiplicities and processing weights from input specified files 
(also listed). Input data sets can be concatenated for multiple runs. 

The first step in the program is the creation of a scene weighted by 

2 

a (sine) antenna power pattern with the necessary null-to-null FOV, and the 
Jacobian factor. This is stored in array SCENE corresponding to K n in 
section V. The unweighted scene has a constant nominal brightness 
temperature. Receiver responses are then calculated according to section V 
and stored in complex array R. The weightings W nm are based on an ideal 
rectangular bandpass filter. A scene estimate SHAT is generated by Fourier 
transforming R. FOURT is a Fourier transform subroutine from the ORI program 
library. The difference between SHAT and SCENE and the fractional difference 
are computed as measures of estimation error. The fractional difference 



TABLE VI - 1 

APSYNSIM Input Variab'es 


Name 

Meaning 

Format 

RIDNAM 

Run ID name 

A20 

N 

Number of baselines 

13 


(power of 2) 


SWDTH 

Scene full width (degrees) 

F8.2 

FOV 

Null-to-null FOV (degrees) 

F8.2 

FO 

Center frequency (GHz) 

F8.2 

BW 

Bandwidth (MHz) 

F8.2 

TEMP 

Constant scene temperature (K) 

F8.2 

MULTFI 

Multiplicity file name 

A20 

WGTFIL 

Weighting file name 

A20 

OUTFIL 

Output file name 

A20 



{ SCENE-SHAT) /SCENE is a measure of the ranometric error induced since SCENE 
and SHAT both are brightness distrioucion" weighted by the same antenra and 
Jacobian factors. Since the scene nomina temperature is specified, 
multiplying the fractional difference by t-.e noiinal temperature gives the 
radiometric error. 

Output at the various stages of the calculation are displayed on 
printer plots generated by PLOT. In addition, an output file is created with 
the input data as a header and four columns containing the abscissa (scene 
coordinate x) and three ordinates-the test scene, scene and estimate 
difference, and radiometric error. The output file is in a format suitable 
for importation to LOTUS from which graphics can be readily generated. 
Simulation results are discussed in the next section. 

A second main program, named BASELINE, is included here. It reads 
from a keyboard specified input file and writes to a keyboard specified output 
file. The input file contains integer antenna locations in a one-dimensional 
array. The output file lists the multiplicity for each baseline in the array. 
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W=E X p ( 1 5 1 GN+2+P I +-SDRT ( - 1 ) +M / ( 4+MM^ ■. ) ) . C HELD FDR W= IS I GW- SOFT ( - ; 
AMD REPEAT FOR W= I SIGN ‘SORT ( — 1 ) +CC.M JUGATE • 1*1) . 


N0M2T=N0N2+N0N2 
IF AR=N7 WO/NPi 

I F ( I PAR-2 ) 350 , 330 , 320 &tiG II 

IPAR=IPAR/4 OF pf 

GO TO 310 

DO 340 11=1,11 RMG , 2 
DO 340 J 3= 1 1 , N0N2 , NP 1 
DO 340 K1=J3,NT01' ,N0N2T 
h 2=1 1+N0N2 
TEMPR=DATA 0,2) 

TEMPI=DATA (12*1) 

DATA <F2) =DA7 A (1 1 ) -TEMPft 
DATA 0 2+l)=DATAO 1 + 1) -TEMPI 
DATA O l)=DATAO 1 ) +TEMPR 
DATA 0, 1*1 ) =DATA 0,1 + 1 ) +TEMPI 
MMAX=N0N2 

I F ( MMAX -NP2HF ) 370 , 600 , 600 
LMAX=MAXO (N0N2T , MMAX/2) 

I F ( MMAX-N0N2 ) 405 , 405 , 380 

THET A=-TWOP I + FLO AT (NGN 2 ) /FLOAT (4-MMAX) 

I F ( I S I ON ) 400 , 390 , 390 

THETA=-THETA 

WR=COS (THETA) 

WI=SIN (THETA) 

WSTPR=— 2.*WI+WI 
WSTPI=2. *Wft+WI 
DO 570 L=N0N2 , LMAX , N0N2T 
M=L 

IF ( MMAX— NGN2 ) 420 , 420 ,410 

W2R=WR*WR-WI*WI 

W2 I =2 . *WR+WI 

W3F:=W2R*WR-W2 1 *W I 

W3 I =W2R* W I +W2 I +WR 

DO 530 11-1,11 RNG , 2 

DO 530 J3= 1 1 , N0N2 , NP 1 

1 M I N= J 3+ 1 PAR+M 

IF (MMAX-N0N2) 430,430,440 

l'MIN=J3 

1 D I F= I PAR* MM AX 

KSTEF-4+PDIF 

DO 520 K1 =LMI N , NTOT , t 5TE P 
K2=K 1 H DIF 
K3=h 2*) DIF 
K4=F3+FDIF 

I F (MMAX-N0N2 ) 460 , 460 , 480 
U1R-DATA <kl) *DATA(I'2) 

UlI=DATAO 1 + 1 ) *DAT A ( K2+ 1 ) 

U2R=DATAO,3) +DhFA(K4) 

U2I=DATA (K3+1 ) +DATA (I 4+1) 

U3R=DATAO 1)-DATA(K2) 

U3 1 =DAT A ( ( 1 + 1) -DATA 0,2+1) 

I ■- ( I S I GN ) 4 70 ,475, 475 
1_..jR=DATA 0,3+1 ) -DATA (1,4+1 > 


Sfe MGE » 

POOR quality 


Microsoft FC)f\'[ fiAt J 


5 


„j 

5 

5 

5 

5 

5 


Li nett 

169 

170 

171 

172 

173 

174 

175 

176 

177 

178 

179 

180 
181 
182 
183 


184 


1 


475 


480 


5 

185 


5 

186 


5 

187 

490 

rr t 

188 


5 

189 


tr ( 

190 

500 

5 

191 


5 

192 

5 1 0 

5 

193 


5 

194 


5 

195 


5 

196 


5 

197 


b 

198 


5 

199 

520 

4 

200 


4 

201 


4 

202 


4 

20 3 

530 

2 

204 



205 


9 

206 

540 


207 


2 

208 


2 

209 


2 

210 

550 

O 

211 


n 

212 


2 

213 

560 


214 

565 

'“i 

215 


n 

216 

570 

l 

217 


l 

2)8 


l 

219 


l 

220 

C 

l 

221 

C 

l 

»-» 

C 

i 

-f 

C 

i 

224 

r 


7 

U4 I=DATA (1 4) -DATA (I 3) 

GO TO 510 

U4R=DAT A 0 4+ i ) —DATA <) 3+1) 

U4 I -DA FA (f 3) -DAI AO 
GO TO 510 
T2R-W2R+DATA (l 2) -W2I*DATA (K2+1 ) 

T2 1 =W2R*DAT A ( t'2+ 1 ) +W2I+DATA <) 2) 
T3R=WR*DATA<L3> -W I* DATA (K3+ 1 > 

T3I=WR+DATA (L3+1 ) +WI+DATA U 3) 
T4R=W3R*DATA'K4) -W3I+DATAO 4+1 ) 
T4I=W3R+DATA O 4+ 1 ) +W3 I +DATA O 4) 

U1R=DATA (hi ) +T2R 
UlI=DATA(kl+l)+T2I 
U2R— T3R+T 4R 
U2 1 =T3 1 +T 4 I 
U3R=DATAU D-T2R 
U3I =DAT A (LI + 1 ) — T2 1 
I F ( I S I GN ) 490 , 500 , 500 
U4R=T3 1 — T 4 1 
U4 I =T 4R-T3R 
GO TO 510 
U4R=T4T-T3l 
U4I=T3R-T4R 
DATA (1 1)=U1R+U2R 
DATA (Ll+1 ) =U1 I+U2I 
DATA (F2) =U3R+U4R 
DAT A ( L2+ 1 > =U3 1 +U4 1 
DATA (H'3) =U1R— U2R 
DATA O 3+ 1 ) =U1 I— U2I 
DATA (K4) =U3F:-U4R 
DAT A ( L4+ 1 ) =U3 1 — U4 1 
KMIN=4* (KMIN-J3) +J3 
1 D I F=1 STEP 

I F a D IF-NF‘2 ) 450 , 530 , 530 

CONTINUE 

M=MMAX-M 

I F ( I S I GN ) 540 , 550 , 550 

TEMPR=WR 

WR=-W I 

WI=-TEMPR 

GO TO 56 u 

TEMPR=WR 

WR=WI 

W I =TEMPR 

I F (M-LMAX ) 565, 565 ,410 
TEMPR=WR 

WR=WR*WSTPR-W I * WSTP I +WR 
l*J I -W I +WSTPR+TEHPR+WSTF' 1 +W I 
I F'AR=3— I PAR 
MMAX =MMAX +MMAX 
GO TO 360 

MAIN LOOP FOR FACTORS NOT EQUAL TO TWO . 
W-EXP ( ISIGN*2-*PI~SQRT <- 1 > * < J2-1 ) - ( J 1 -J2) / 
PERFORM A FOURIER TRANSFORM OF LENGTH IFf- 
CONJ UGATE SYMMETR 1 E 3 . 
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ORIGINAL PAGE IS 
OF POOR QUALITY 


APPLY THE TWIDDLE FACTOR 
(NP2*IFP1>>, THEN 
CT (IF), MAILING USE OF 


| 3 S 


f ■ c> a tc 


— C *■ 


D 

Lint- ** 

1 

7 Mi eras:- a •' 

1 

— . '-VCT 

C 


J 

226 

d>00 

IF (NTW0-NP2) 603 , 700 , 700 

f 

?T1 

605 

IFP1=N0N2 

1 

228 


I F= 1 

1 

229 


NP 1HF =NP 1/2 

1 

230 

6 1 0 

I FP2= I FP 1 / I FACT ( IF) 

1 

231 


J 1RMG-MP2 

1 

232 


IF < I CASE-3) 612,61 1 ,612 

1 

933 

61 1 

J 1RNG= (NP2+IFP1) /2 

1 

234 


J2STF-NP2/IFACT ( IF) 

1 

23 5 


J 1RG2= (J2STFM-IFP2) 72 

1 

236 

612 

J2HIN=1+IFP2 

1 

237 


IF (IFP1-NP2) 615,640,640 

1 

238 

615 

DO 635 J2=J2MIN, IFP1 , IFP2 

2 

239 


THETA=-TW0PI -*FL0AT ( J2-1 > /FLOAT (NP2) 

2 

240 


I F ( I S I BN) 625 , 620 , 620 

n 

241 

620 

THETA- -THETA 

2 

242 

625 

S I NTH=S IN < THETA/2. ) 

n 

243 


WSTPR=-2.-*SINTH*SINTH 


244 


WSTP 1=0 IN (THETA) 


243 


WR=W5TPR+1 . 

'—i 

246 


V>J I =WSTP I 


247 


J 1MIN=J2+IFF’1 

n 

240 


DO 635 J 1 = J 1 M I N , J 1 RNG , I FP 1 

-r 

249 


I 1 MAX=J 1 + I 1RNG-2 

3 

250 


DO 630 1 1 = J 1 , UMAX, 2 

4 

251 


DO 630 13=11 ,NT0T,NP2 

5 

■i— 


J3MAX= I 3+ I FP2— NP 1 

k 5 

253 


DO 630 J3= I 3 , J 3MAX , NP 1 


254 


TEMPR-DATA ( J3) 

6 

255 


DATA (J3) =DATA ( J3) -tWR-DATA ( J3+1 ) *l‘JI 

6 

256 

630 

DATA ( J3+1 ) =TEMPR*W I + DAT A ( J3+ 1 ) *WR 

3 

257 


TEMPR-WR 

_> 

258 


WR=WR* WSTPR-W I *WSTP I +WF: 

T, 

259 

635 

W I =TEMPR itWSTP I +W I ^WSTPR+W I 

1 

260 

640 

THETA=-TUJOPI /FLOAT ( I FACT ( IF) ) 

1 

261 


IF (ISIGN) 650,645, 645 

1 

262 

645 

THETA=-THETA 

1 

263 

650 

SINTH=SIN (THETA/2. ) 

1 

264 


WSTPR=-2. *SINTH*SINTH 

1 

265 


WSTP I “SIN (THETA) 

1 

266 


LSTEF'=2*N/ IFACT (IF) 

l 

267 


KRANG=I STEP* (IFACT (IF) /2>+l 

1 

268 


DO 698 11=1,1 1RNG , 2 

n 

269 


DO 698 13=1 1 ,NT0T,!MP2 

3 

270 


DO 690 LMIN=1 ,) RANG,) ETEP 

4 

271 


J 1MAX=I3+J 1RNG-IFP1 

4 

272 


DO 680 J1=I7,J1 MAX , I FP 1 

IT 

»_» 

273 


J3MAX=J 1 +• IFP2-NF' 1 

C7 

274 


DO 680 J3=J 1 , J3MAX , NF' 1 

6 

275 


J 2MAX=J3f IFF' 1 — IFF' 2 

6 

27o 


K=KMIN+ < J3-J1+ (J1-I3) /IFACT < 1 ^) ) /NP1HF 

6 

277 


IF (1 MIN-1 ) 655 , 655 ,665 

■•D 

270 

655 

SUMR-O. 


279 


SUM 1=0. 

i 43 

280 


DO 660 J2=J3 , J2MAX , IFF’ 2 



/ 



D 

Li no h 

i 

i 

7 

281 


F _ • 1 F:=;‘ UMR+ DATA U2) 

7 

282 

6 c 0 

ELM I — LUf-l I h-DAT A (J2+1) 


203 


a > =sut r 

T h 

284 


o / 

ii 

/*. 

4 

6 

285 


L-. TO r 80 

6 

286 

665 

1 CGNJ =~ +2* IN-KMIN+1) 

6 

287 


J 2= J 2MAX 

6 

288 


SUMR=DATA ( J2) 

6 

289 


SUM I =DAT A < J 2+ 1 ) 

6 

290 


C"-.DSP:=0 . 

6 

291 


OLDS 1=0. 

6 

292 


J 2= J 2- 1 FP2 

6 

293 

670 

TEMPR=SUMR 

6 

294 


TEMP I =SUMI 

6 

295 


SUMR=TWOWR*SUMR-OLDSR+DATA ( J 2 > 

6 

296 


SUM I =T WOWR^SUM I —OLDS I +-DA T A (J2+ 

6 

297 


OLDSR=TEMPR 

6 

29S 


OLDS I =TEMP I 

6 

299 


J2=J2-IFP2 

6 

300 


IF (J2-J3 >675,675, 670 


llicraec!) t FCiF.7ftA.N77 


r ut_ 

0 2 -• 

>0 Mi-rcn 1 


TEMPR=WR*SUMF:-OLDSR+DAT A ( J2 ) 

TE M F‘ I=WI ♦SUM I 

WQR1 (^ ) =TEMPR- 7 EHP I 

W0R1 ( h CON J ) — TEMPR+T EMF I 

TEMPR=WR+SUMI-OLDSI+DATA (J2+1 > 

TEMP I = W I *SUMR 
WOR1 (ft+l)=TEMPR+TEMF'I 
W0R1 (KCONJ+1 ) =TEMPR— TEMPI 
CONTINUE 

IF (KMIN-l ) 685,605,686 

WR=WSTPR+1. ' . 

WI=WSTPI 
GO TO 690 
TEMPR-WR 

WR=WR*WSTPR-W I *WSTP I +WR 
W I =TEMPR*WSTP I +W I *WSTPR+W I 
TWO W R = W R + W R 

IF (ICASE-3 >692,691 ,692 

IF ( I FP1-NP2) 695,692,692 

ft=l • 

I 2MAX = I 3+NP2—NP 1 
DO 693 12=13, I2MAX ,NP1 
DATA (12) =W0R1 (K> 

DATA! 12+1 )=WORF U +1) 

K=fs+2 

GO TO 698 

COMPLETE A REAL TRANSFORM IN THE 1ST DIMENSION, N ODD, BY CON- 
JUGATE SYMMETRIES AT EACH STAGE. 

J 7MAX= I 3+ I FP2-NP 1 

DO 697 J 3= 1 3 , J 3MAX , NP 1 

J2MAX=J3+NP2-J2STP 

DO 697 J2=J3, J2MAX , J2STP 

J I MAX=J2+J i RG2- IFP2 

: i CNJ=J3+ J2MAX +J2ST P-J2 





Lx nett 1 
337 
339 

339 

340 

341 

342 696 

343 

344 697 

345 698 

346 

347 

348 

349 C 

350 C 

351 C 

352 C 

35 3 700 

354 701 



358 702 

359 703 

360 

36 1 

362 

363 

364 

365 

366 

367 710 

368 

369 

370 

371 

372 

373 

374 

375 

376 

377 

378 

379 720 

380 

381 

382 

383 

384 

385 725 

386 734) 

387 731 

388 735 

389 740 
394) 

391 

392 


F'aqi? 
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CO c-n J 1=J2. J1MAX , IFP2 

' = i + : i - 1 3 

0-7 A Jl)=WOR» (3 > 

~m7h J i + 1 > = ! 0CJP3 (3+1 ) 

IF x J x-J2) 697 ,c97,696 

data (jicmj) =woRf <30 ORIGINAL PAGE IS 

DATA x J 1CNJ+1 ) =-W0R3 (3+1) OF POOR QUALITY 

J 1 CN J — J 1 CN J — I FF'2 

CONTINUE 

IF= IF+ 1 

IFP1 = IFF'2 

I F ( I F P 1 -NP 1 ) 700 , 700 ,610 

COMPLETE A REAL TRANSFORM IN THE 1ST DIMENSION, N EVEN, BY CON- 
JUGATE SYMMETRIES. 

GO TO ( 900 , 000 , 900 , 70 1 ) , I CASE 

NHALF=N 

N=N+N 

THET A--TWOP I /FLOAT (N) 

I F ( I S I GN ) 703 , 702 , 702 
THETA- -7 HE 'F A 
SINTH=SIN (THETA/2. ) 

WSTPR=-2. V S I WTH*S INTH 
WSTPI=S IN (THETA) 

WR=WSTPR*1 . 

W I =WSTP I 
I M I N— 3 

JMIN=2*NHALF-1 
GO TO 725 
J=JMI N 

DO 720 1= IMIN , NTOT , NP2 
S'JMR= (DATA ( I ) +DATA ( J ) ) /2. 

SUMI = (DATA (I + D+DATAXJ + l) ) /2. 

DI FR= (DATA ( I ) — DATA ( J ) )/2. 

DIFI= (DATA ( 1+1 ) -DATA <J + 1> ) /2. 

TEMPR=WR*SUM I +W I *DI FR 
TEMP I =W I -x-SUMI -WR*D I FR 
DATA ( I ) =SUMR+TEMPR 
DATA < I + 1 ) =D I F I +TEMP I 
DATA ( J ) =SUMR-TEMPR 
DAT A ( J + 1 ) =— D I F I +TEMP I 
J= J +NP2 
IMIN=IMIN+2 
JMIN=JMIN-2 
TEMPF;— WF: 

WR=WR x WSTPR-W I * WST P I -t WR 
W I =TEMPR k WST P I -*-W I - WSTF P+l*J I 
I F ( I MIN- JM I N ) 7 1 0 , 730 . 74 0 
IF ( IS IGN> 731 ,740, 7*i4> 

DO 735 I = I M I N , NTOT , NP2 
DATA ( I + 1 ) =— DA TA(IH) 

NP2=NP2+NP2 
NTOT =NTOT -* NTOT 
J=NT0T+1 
1 MAX —34 TOT / 2-*- 1 





1 0 L 

lie ri 

X 

7 H c 

» 1 

39 7. 

745 

I M I N= 1 1-1 A X —2* NHALF 

1 1 

3^4 


I - I M I N 

m 

7^5 


GO TO 755 


3 C> 6 

750 

D*TA ( J) = DATA ( I ) 

1 

397 


D-TA <J + 1 ) =-DAT A ( 1 + 1 ) 

1 

39 8 

755 

1=1+2 

1 

399 


J = J -2 

1 

400 


I F ( I - 1 MA X ) 750 , 760 , 760 

1 

40 1 

760 

DATA ( J ) =DATA ( I M I N ) -DATA ( I M I N+ 1 ) 

1 

402 


DAT A ( J+ 1 ) =0. 

1 

*03 


IF (I-J) 770, 780, 780 

1 

404 

765 

DATA ( J ) =DATA < I ) 

1 

405 


DAT A ( J + 1 ) =DAT A < I + 1 ) 

1 

406 

770 

1=1-2 

i i 

407 


J = J — 2 

i 

400 


IF ( I -I MIN) 775,775, 765 

i 

409 

775 

DATA ( J ) =DAT A ( IMI N) +DAT A (IMIN+1 ) 

i 

410 


DATA < J + l ) =0. 

i 

411 


INAX=IMIN 

i 

412 


GG TO 745 

i 

413 

780 

DATA < 1 > =DAT A ( 1 > +DATA (2) 

i 

414 


D^-TA (2) =0. 

i 

415 


GG TO 900 

i 

416 

C 


i 

417 

C 

COMPLETE A REAL TRANSFORM FOR THE 

i 

410 

c 

CONJUGATE SYMMETRIES. 

i 

419 

c 


i 

420 

800 

I F ( I 1 RNG-NP 1 ) 805 , 900 , 900 

^ i 

421 

805 

DO 860 13=1 ,NTOT ,NP2 

m 2 

422 


1 2M A X = 1 3+NP2— NF' 1 

W -2 

423 


DO 860 12=13, I2MAX,NP1 

3 

424 


IMI N= I 2+ I 1 RNG 

_> 

425 


IMAX=I2+NPl-2 

3 

426 


JMAX=2* 1 3+NP 1 — I M I N 

T 

427 


IF (12-13)820,820,810 


428 

810 

JMAX=JMAX+NF2 

3 

429 

820 

IF ( IDIM-2) 850,850,830 

3 

430 

830 

J=JMAX+NP0 

*T ( 

431 


DO 840 I = I MI N , IMAX , 2 

4 

432 


DAT A ( I ) =DATA ( J ) 

4 

433 


DATA ( 1 + 1 ) =— DATA ( J + l ) 

4 

434 

840 

J = J -2 

3 

435 

850 

J = JMAX 

3 

436 


DO 860 I=IMIN, IMAX ,NPO 

4 

437 


DATA ( 1 ) =DATA < J ) 

4 

438 


DATA ( I + 1 ) =-DATA ( J + 1 ) 

4 

439 

860 

J =J— NF'O 

i 

440 

C 


1 

441 

C 

END OF LOOP ON EACH DIMENSION 

1 

442 

C 


X 1 

443 

900 

N-'0=NP1 

9 1 

44 4 


N C 1=NP2 

I i 

445 

910 

N- REV=N 

1 

446 

920 

RETURN 

F 

^,17 


END 
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IS\ 



I Line# 1 


/ 



C R 
DATA 

d:fi 

DIFR 

FLOAT 

I 

I I 

I lfiAX 
I1RNG 
12 

I2MAX 

13 

I CASE 
IDIM 
IDIV 
IF 

IF ACT 

I FORM 

IFF 1 

I FP2 

IMAX 

IMIN 

I PAR 

IQUOT 

IF'EM 

ISIGN 

J 

J1 

J1CNJ 
J 1 MAX 
J1MIN 
J1RG2 
J 1 RNG 
J2 

J2MAX 

J2MIN 

J2STP 

J3 

J3MAX 

JMAX 

J M I N 

k 

kl 

K2 

K3 

K4 

1 COM3 
i DIF 
KM IN 
KRANG 
KSTEP 

i_ 

L.MAX 
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Type 

Of * set 

REAL 

0 

REAL 

704 

REAL 

700 

INTEGERS 

208 

INTEGER-4 

252 

INTEGER-4 

248 

INTEGER-4 

220 

I NTE6ER-4 

236 

INTEGER-4 

640 

INTEGERS 

260 

INTEGER-4 

204 

.•INTEGER- 4 

152 

INTEGER* 4 

188 

INTEGER* 4 

18^ 

INTEGER-4 

16 

INTEGER- 4 

16 

INTEGER - A 

472 

INTEGER-4 

480 

INTEGER-4 

716 

INTEGER-4 

684 

I NT EGER- 4 

288 

INTEGER-4 

192 

I NT EGER- 4 

196 

INTEGER-4 

12 

INTEGER-4 

212 

INTEGER-4 

520 

INTEGER- 4 

66G 

INTEGER-4 

580 

INTEGER-4 

516 

INTEGER-4 

492 

INTEGER-4 

484 

INTEGER-4 

5 0 0 

INTEGER-4 

600 

INTEGER-4 

496 

INTEGER-4 

488 

INTEGER-4 

°72 

INTEGER-4 

544 

INTEGER-4 

736 

INTEGER-4 

688 

INTEGER-4 

604 

INTEGER-4 

304 

INTEGER-4 

316 

INTEGER- 4 

408 

INTEGER- 4 

412 

INTEGER- 4 

624 

INTEGER- 4 

392 

INTEGER-4 

388 

INTEGER* 4 

556 

INTEGER- 4 

386 

INTEGER- 4 

348 

IhPEGEP— 4 

3~j<i 


Class 

INTRINSIC 


INTRINSIC 


\to 


i 



0 Lini-t; 2 


1-1 

MAXO 

INTEGER-M 

176 

w* 

I NT EGER* 4 

320 

INTEGER- 4 

168 

■JOIN 

INTEGERS 

8 

NHALF 

INTEGER* 4 

660 

NN 

INI EGER*4 

4 

NON 2 

INIEGER+4 

200 

N0N2T 

INTEGER+4 

284 

NPO 

INIEGER*4 

224 

NP1 

INTEGERS 

160 

NP 1 HP 

INTE6ER*4 

476 

NP2 

INTEGERS 

172 

NP2HF 

INTEGER* 4 

232 

NPREV 

INTEGER*4 

228 

NTOT 

INTEGER*4 

148 

NTWO 

INTEGER*4 

180 

OLDSI 

REAL 

632 

OL.DSR 

SIN 

REAL 

628 

SI NTH 

REAL 

512 

SLIM I 

REAL 

612 

SUMP 

REAL 

608 

T2I 

REAL 

452 

T2R 

REAL 

448 

T3I 

REAL 

460 

T3R 

REAL 

456 

T4I 

REAL 

468 

TAR 

REAL 

464 

Itempi 

REAL 

280 

TEMFR 

REAL 

276 

THETA 

REAL 

328 

TWOP I 

REAL 

144 

TWO NR 

REAL 

636 

U1 1 

REAL 

420 

U1R 

REAL 

416 

U2I 

REAL 

428 

U2R 

REAL 

424 

U3I 

REAL 

436 

U3R 

REAL 

432 

U4I 

REAL 

444 

U4R 

REAL 

440 

W2I 

REAL 

364 

W2R 

REAL 

360 

l*J3 1 

REAL 

372 

W3R 

REAL 

368 

WI 

REAL 

336 

WORK 

REAL 

20 

WR 

REAL 

332 

W5TPI 

REAL 

344 

WSTFR 

REAL 

340 

Name 

T ype 

Sice 


INTRINSIC 


INTRINSIC 


Cl ass 


Hi c i o 
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VII. SIMULATION RESULTS 


In this section we present our simulation results using a simple but 
realistic sensor as an example. The sensor parameters are listed in table 
VII-1. A scene width of 90° and a nominal temperature of 250K are used. 


The sensor resolution is approximately the scene sampling interval 
(90°/32). At an altitude of 300 km this corresponds to 15 Km resolution 
with an antenna array dimension given approximately by 



C . 1 ( 2.CG k 5 

I 5 





The value .21 m is the wavelength. There are no grating lobes as the value 





(. 2 1 /VI ) 

15 
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Figure VII-1 shows the constant test scene weighted by the antenna 
power and Jacobian factors as a function of scene coordinate (sineof across 
track angle). Figures VI 1-2 to 7 show the weighted errors in the scene 
estimates. By this we mean the error before the antenna and Jacobian factor 
are removed. Figure 2 is for zero bandwidth, in this case the Fourier 
processing is perfect and only random round-off error is present. Figures VII 
- 3 to 7 have bandwidths of 10, 20, 30, 40, and 50 MHz respectively. As the 
bandwidth increases the error amplitude envelope increases. 



TABLE VII - 1 
Sensor Parameters 


Name 


Parameter 


Value 


N 

Number of baselines 

FOV 

Field-of-view 

FO 

Center frequency 

BW 

Bandwidth 


Uniform multiplicities 
Uniform processing weights 


32 

60° 

1.4 GHz 

0, 10, 20, 30, 
40, 50 MHz 


150 



The most important information is contained in figures VII - 8 to 12 
which show the radiometric (or unweighted) errors. These errors are small but 
nonzero near nadir (scene coordinate 0). Off -nadir they become large, 
especially near the pattern nulls. In general, a region around nadir exists 
where radiometric errors are uniformly tolerable. Defining this region to 
have radiometric errors less than say 1.0 K specifies the useful field-of-view 
of this sensor. In figure VII - 13 we have plotted useful field-of-view 
versus bandwidth by interpolating the scene coordinate for which the 
radiometric error exceeds 1.0 K. At zero bandwidth the entire 6(P 
field-of-view is usable, but this diminishes as the bandwidth increases. 
Finally, at 50 MHz there is no usable field-of-view since nadir radiometric 
errors exceed 1.0 K. The bandwidth of 30 MHz (just greater than that 
allocated for passive use) has a useful field-of-view of about 3(JP by this 
definition. Obviously there is a tradeoff here between radiometric errors 
induced by this effect and the receiver radiometric sensitivity, since the 
former increases with bandwidth and the latter decreases. 
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Figure VII-3. Scene Estimate Error, BW = 10 MHz. 
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Figure VI 1-4. Scene Estimate Error, BU * 20 MHz. 
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Figure VII -6 


Scene Estimate Error, BW = 40 MHz 
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Figure VII-7. Scene Estimate Error, BW = 50 MHz. 
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Figure VI 1-8. Radiometric Error, BW = 10 MHz 
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Figure VI 1-9. Radiometric Error, BW = 20 MHz 
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Figure V I I — 1 0« Radiometric Error, 


BW = 30 MHz. 
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Figure VII-11. Radiometric Error, BW = 40 MHz. 
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Figure VII-12. Radiometric Error, BW = 50 MHz. 
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Figure VII-13. Useful Field-of-View. 
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V. CONCLUSIONS AND AREAS FOR FURTHER STUDY 


Aperture synthesis may be an attractive approach to perform the soil 
moisture measurement mission. It can have both a resolution and a sensitivity 
advantage over real aperture systems. Further study is needed especially in 
the areas of: 
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1) Technologies 

a) antenna design 

b) signal distribution 

c) digital correlator design 

2) Calibration 

3) Antenna and array structural tolerances 

4) Trade-offs between pure and hybrid aperture synthesis techniques 

5) Trade-offs between real and synthetic aperture systems (i.e., 
when does the cost of electronics exceed the savings from 
antenna area). 

6) Failure modes and reconfigurability 

7) Adapting of the spacecraft concepts to aircraft prototypes. 
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